Function Definition:

R = 8.314
if CalUnits = true then R = R / 4.18

T = Temp(i) + 273.15

Tm0 = Tm + 273.15
dH0 = dHi * 1000
K = exp(-dH0*(1-T/Tm0)/(R*T))
Cp=(dH0 * dH0)/(R*T*T)*exp(dH0*(-Tm0+T)/(Tm0*T*R))/((1+K)*(1+K))*aw
Y = Cp / 1000
i = i + 1



Globals:

dim i
i = 0